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Abstract 

We consider the behavior of the overlap of m[> 2) paths at the spin glass 
transition for a directed polymer in a random medium. We show that an 
infinite number of exponents is required to describe these overlaps. This is 
done in an e = d — 2 expansion without using the replica trick. 
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When disorder induces a new thermodynamic phase not found in the pure system, the 
description of the transition itself poses new challenges. This is exemplified by the efforts 
made to understand the spin glass transition [00. Other examples are magnets in random 
fields 0, polymers in disordered media f^-j^ or with random interactions p|, etc. Several 
concepts emerged from the solution of the infinite range "mean field" spin glass problem using 
the replica trick, but their validity for finite dimensional systems, remains to be established 

Em- 

The most important concept in the replica approach is that of the overlap which plays the 
role of the order parameter for the transition The overlap purports to characterize the 
rugged free energy landscape in the spin glass phase through the statistics of the pairwise 
common configurations in the various minima. In principle, a many valley free energy 
surface would require higher order overlaps for a more detailed description [p^p!2| . These 
are the overlap of say, three or more (m) minima, to be called m- overlap" . All of these 
overlaps vanish at the transition point as ^\ — T |^'" for T — > T^— ,because the free 
energy surface goes over to a smooth one in the high temperature phase. Now, how many 
exponents are needed to describe these overlaps? The answer is one for infinite range or 
infinite dimensional models ll|Jll|, justifying the use of the pair overlap, q2, as the sole order 
parameter in the replica approach (with j32 as the order parameter exponent) |TB|p!^ . What 
about finite dimensional systems? The question assumes importance because the number 
of independent exponents tells us the number of quantities one requires to characterize the 
transition. Alas, so little is known for the spin glass problem 0. 

In this paper, we develop a method to calculate these m-overlaps for arbitrary m in an 
e = d — 2 expansion for the spin glass transition of a directed polymer in a random medium 
(DPRM) without using the replica trick. A d + 1 dimensional DP is a string stretched 
in a preferred direction and with free fluctuations in the transverse d dimensions. In a 
random medium, the gain in the potential energy from randomness can win over the "random 
walk" entropy, producing a disorder dominated 'super diffusive" phase [p|-[7| JT^JTB[ | . For 
d > 2, there is a transition from a low temperature, strong disorder, spin glass type phase 
to a pure type phase The transition is described by an unstable fixed point 

(~ 0(e)) in a renormalization group (RG) approach via the mapping [|1^] to a nonlinear 
noisy stochastic equation for the free energy (Kardar-Parisi-Zhang equation) PJ^. The 
simplicity of the model, nontrivial solutions in many cases [pO| , p!T] , p!2| , and the possibility of 
studying various fundamental questions related to disorder systems in general, make 

the strong disorder phase and the transition a topic of paramount importance In fact, 
various techniques have been used for this purpose 



as. 



e.g., Bethe ansatz for = 1 pO 



dynamic renormalization group scaling theory |TB], numerical simulations [23|, 

mode coupling p6[, 1/d expansion p^, on a Cayley tree |[Tl|], on hierarchical lattices [|12| etc. 
Replica symmetry breaking has been tried by some and vehemently opposed 

by a few others [pO| , pl[1 . 

For DP, the overlap (m-overlap) describes the fraction of paths common to two (m) 
minimum (free) energy paths, and is therefore equal to the fraction of the paths two (m) 

. The two 



polymers go together when placed in the same random medium [p!5| , p!q , p!l| , pr^ , P^ ,[ 
(m) polymers act as the replicas of the one chain system. In the spin glass phase, if there is 



only one minimum free energy path, then both the chains would follow the same path, an 
attraction induced by the disorder. In case there are many valleys, then the chains can get 
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separated by hopping to a neighboring valley. Thus overlaps contain information about the 
valley structure. 

The m-overlaps, g^, have been calculated for a DPRM on a Cayley tree and on hierar- 
chical lattices [0,|ll|] • The Cayley tree problem can be thought of as an infinite dimensional 
case while the hierarchical lattices are definitely finite dimensional with tunable dimension- 
ality. For the Cayley tree problem, closed form expressions for show that (3m = 1 for 
all m [|TT|. For hierarchical lattices, there is a critical dimension above which a transition 
takes place |jl2|]. Numerically the transition temperature has been obtained by locating the 



temperature where q2 and vanish. Nothing, unfortunately, is known for (3m- A linear 
dependence of (3m on m, in a multifractal analysis, would also mean that only one exponent 
is needed, as e.g. for pure noninteracting Gaussian chains, (see below). The crucial question 
is, therefore, whether such linearity is maintained for finite dimensional systems. The answer 
we find is no. 

In this paper, we use the continuum approach. The polymer is described by the Hamil- 
tonian 

'7 _.2/ ^ , ^ 



H = dr 







2 x^(r) + -nx(r),r) 



(1) 



where x(t) is the d dimensional transverse spatial coordinate of the DP at the contour 
length t, and x(r) = (ix(r) / dr. The first term on the right hand side represents the entropic 
fluctuations of a free Gaussian chain with 7 as the bare line tension. V corresponds to a space 
and time dependent random potential seen by the chain, and the amplitude A/27 is chosen 
for convenience. The random potential is taken to be uncorrelated, normally distributed 
with zero mean and V {x, t)V {y , t') = 2A 5(x — y)5(r — r') where the overbar stands for 
disorder averaging. 

A formal way to define the m-overlap is to put m chains in the system and take 



I rt /m-l 



V)) , (2) 



where Xj j+i(r) = Xj(r) — Xj+i(r), and (...) stands for thermal average for a realization. One 
way of computing is to couple the m chains or replicas with a weak m-body interaction. 
Then overlap would follow from an appropriate derivative of the total free energy with 
respect to the hypothetical coupling constant (see below). This procedure was adopted for 
the overlap (^2) in the numerical work of Mezard in 1+1 dimensions and by one of us in 
a 1-loop RG approach (see also Ref. |]52|). We generalize the method of Ref for qm- 
The RG analysis is geared towards calculating the scaling exponent of the coupling constant. 
A judicious use of finite size scaling, as explained below, then gives us the exponent /3m- 

With the definition of the m-overlap in Eq. ^ we consider an m chain interacting 
Hamiltonian 

m „i m— 1 

nm = Y.H^ + iX/2^)vm / dr [] 5(x,,+i(r)). (3) 
t=l -^o i=l 

where Hi is the Hamiltonian of Eq. |I| for the ith polymer. Defining the quenched free energy 
fmi^myt) = In wherc Zm is the partition function for Hm, the m-overlap is obtained as 

Qm = ^ dfm{Vm^t) /dVm\vm=0- 
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Our interest is in the scaling part of the free energy, ^ t>^/^f(i;^t-<^™/^), where (pm 
is the scahng exponent for Vm, and x and z are the single chain free energy fluctuation 
(A/ ~ t^/^) and dynamic (size x ~ t^^^) exponents ||2^j5[]. We have verified that, as 



expected, there is no change in the single chain exponent. Taking derivative then gives 
Im ~ t^"", with = {x ~ (l^m — z)/ z. Incidentally, the system is taken to be infinite in 
extent in all the d transverse directions and is of length t in the preferred direction. This 
form of Qm can, therefore, be treated as a finite size scaling form [Q. Now, the transition 
takes place only in the thermodynamic limit of t —>■ oo. In that limit, for d > 2, there is 
a diverging length scale with exponent u, ~| T — \~'^ , parallel to the special t-like 
direction ||5|,p!^. Finite size scaling suggests a scaling form q.^ = t~^"^^'^g{t/^\\)- Therefore, 



right at the critical point (i.e., the unstable fixed point in RG), Qm ~ t"^"/'^. A comparison 
then yields Pm = —t^^m- Remember, that u is strictly m independent. Our strategy is 
therefore to calculate 0m- 

Define /i({xj},t) = (27/A) In Z({xj}, t), where Z({xj},t) is the partition function for 
chains with end points at {x^}, all starting at the origin. This h satisfies the equation 




lV'^h + ^{V,hr^+go, (4) 

where 5^0 = Yl]Li ^(xj, t) +Vm Y[Y=i ^(^j This equation, though resembling a higher 

(md) dimensional KPZ equation [^], is actually not so for the peculiar noise term. We prefer 
this equation to Eqs. |T] or ^ because an RG can be implemented with the nonlinear term 
A as perturbation. This is different from a perturbation in the random potential around 
the Gaussian chains. Moreover, the equation describes the free energy and averaging h will 
naturally give the quenched average free energy without any recourse to the replica trick. 

A scaling of x ^ 6x, t —>■ b^t, then shows that, in the absence of nonlinearity A, Vm 
_^ j^z-{m-i)d~Xy^_ An anomalous part r]m would creep in through renormalization when the 
nonlinearity A is present. The crossover exponent is therefore 0^ = —[z — {Tn — l)d — x + Vm]- 
This gives 

Pm = v[d{m - 1) + r]J^/z. (5) 

The m dependence, apart from the Gaussian one, therefore comes from r]m- At the Gaussian 
level the exponents depend linearly on m. 

The formal solution of Eq. ^, in the Fourier space (kj,u;) conjugate to (xj,t), can be 
written as 

h{{\^j],uj) = Go({k,},a;)(7o - (A/2)Go({k,}, cu) x 

(6) 



^.^ (^E^.j /^({p,},fi)/^({^•-p,},cu-^]), 



where Go({kj},ci;) = [7 ~ iuj\~^ represents the bare m particle propagator (Green's 

function) (see Fig. 1), and go is the Fourier transform of go. A shorthand notation is used, 
viz, Vj = Pj ■ (kj — pj), /|.| = (277)""^°'"-'^/ dVlYYlL^^dpi, and {kj} to represent all the m 
k vectors. To tame possible divergences, we put an upper cutoff A (=1) for all p integrals. 
This cut off actually comes from a short distance cutoff in real space. 
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We now use a dynamic renormalization group approach to determine the behavior of 
Vm for large length scales by integrating out fluctuations on smaller scales. This is based 
on Eq. ^ The procedure is well documented especially for the KPZ equation and for 
the m = 2 case of Eq. ^ |2^. The basic idea is to (i) integrate out the fluctuations at the 



shortest scale by taking a slice Ae~^'' < p < A from the momentum integral, (ii) absorb 
it in the coupling constant, (iii) and rescale all the momenta etc to get back the original 
cutoff A. The resulting changes are then absorbed by renormalizing the parameters of the 
problem. In the limit 5/ — > 0, these changes are expressed in terms of differential equations 
(recursion relations) that tell us the flow of the parameters as we go to longer length scales. 
Special care is needed, for the problem in hand, to keep track of the momenta indices that 
get coupled by the noise. These interchain connections produce the necessary anomalous 
part in the RG equation. We skip the algebraic details. 

The flow of the disorder is described in terms of the dimensionless coupling U = 
KdX^A/{2-f^) where = {27r)-'^Sd, Sd being the surface area of the unit d dimensional 
sphere, and can be found in Ref. 0. We have verified that this single chain equation is 
recovered from Eq. ^ and is independent of m. Let us recapitulate that at = 2, [/ is 
marginally relevant. This leads to a new critical point for d > 2. For d = 2 + e, the unstable 
fixed point U* = 2e corresponds to the spin glass transition point, with u = 2/{d — 2) as the 
length scale exponent ||3^ . 



We concentrate on the renormalization of the coupling Vm- For the long wavelength, 
long time limit, the external wavevectors and frequency are small or zero. In this limit, 
for arbitrary m, the effective coupling constant to one loop order (see Fig. 1) is given by 
(suppressing the zero wave vectors ) 



VmR = V, 




X 



/ p^Go(p,fi)[Go(-p,-f^)]'Go(p,-p,0). (7) 



The recursion relation for Vm follows from Eq. |^ as 

dVrr,. 



dl 



z - X - (m - + 2 j ^]^'»- (8) 



Since our interest is in the crossover exponent for Vm at Vm = 0, this first order (in Vm), 
one loop equation is sufficient for us. Higher loops will generate higher order terms in U 
(and hence e). 

Using the one-loop fixed point value U* = 2e for the transition point, we find rjm = 
—m{m — l)e, which from Eqs. ^ gives 

= Kc[2(m - 1) - (m - Ife + 0{e% (9) 

where Cc = is the size exponent at the transition point. We now see that the linear 
relation between (3^ and m at the Gaussian level is not respected in the first order. In other 
words the exponents f3rn are not linearly dependent on each other, and higher order terms 
are expected to make the interdependence more complicated. Hence the need for an infinite 
number of exponents at the transition point. 
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From the nature of the perturbation series, we see that the change in exponent to first 
order in e is due to the effective two body interaction induced by the disorder. The loop in 
Fig 1 and in Eq. |^ comes when disorder couples two different chains. With 6 correlated noise, 
this means that the two chains are going through the same point in space. The purpose of 
overlap is to count these. It is wellknown that two body interaction changes the reunion 
behavior of m chains, and each m requires a new exponent for reunion p5| , p6| . This is the 
situation here. 

Let us now try to connect this result to a replica analysis. We take n chains in the 
random medium and average the resulting partition function, or equivalently, get the effective 
Hamiltonian for the nth moment of the partition function. The effect of the disorder is to 
couple these chains through a two body interaction The m-overlap, then corresponds to 
the reunion of m chains out of these m, in the limit n — > (the replica trick). At the critical 
point, taking the chains to be Gaussian (since Cc = l/-2c = 1/2) 0,0], the reunion of a 
subset of interacting random walkers can be studied following Ref. . The only difference 



with Ref |3§] is that the relevant fixed point is the unstable one, and the n — > limit can be 
taken, [see, e. g., Ref. ||3^(b)] to get the anomalous part of Eq. ^. This is essentially correct 
but it still needs to be established that the chains are actually Gaussian (not just (c = 1/2). 
These problems are not present in the differential equation approach used in this paper. 

What do all of these mean for the spin glass transition in finite dimensions? One expects 
to to write down (in the n limit) a Landau- Ginzburg type free energy functional with the 
overlaps as the order parameters. There can be two possibilities. (A) One is that the simple 
minded description through the overlap of two copies is not sufficient, and one has to worry 
about the higher overlaps, and in fact an infinite number of them. Even if one starts with the 
Gaussian distribution, renormalization effects will generate the higher overlaps (arbitrary 
distributions generate then in any case). In this situation, the conventional replica approach 
may not be useful. The problem here again may be the interchange of the two limits, 
viz,n and the thermodynamic limit. (B) The other option is that the two copy overlap 
is good enough in the sense that higher m-overlaps are irrelevant. Naively speaking, at the 
transition point with z = 2,x = 0, Vm is irrelevant at ci = 2 for m > 2. However, this does not 
necessarily imply irrelevance of the m-overlaps in the single chain problem. Remember, that 
Vm is a coupling introduced by hand in the many chain Hamiltonian in Eq. ^ to calculate the 
overlaps and does not appear in the description of the single chain problem. This opens up a 
the possibility where under special conditions higher order overlaps can become important as 
in multicritical cases or polymers with higher order composite operators becoming relevant. 
Proper choice of parameters can then lead to multicritical analogs of spin glasses. These 
are not the random version of pure multicritical models but are of inherently different type. 
In both cases, if the m-overlaps are important at the transition, they are expected to be 
so in the spin glass phase also. We hope that this will motivate further detailed numerical 
work to settle this issue. We conclude that, like multifractals, a spin glass transition in finite 
dimensions subsumes an infinite number of independent exponents. This indicates either a 
failure of the simple replica picture in finite dimensions or the possibility of highly complex 
spin glass phases. 

We thank A. Baumgartner for discussions and hospitality at KFA-IFF. The visit of the 
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FIGURES 



FIG. 1. Diagrammatic representation of the parameters (a) and the solution (b). (See also 
Figs 1 and 2 of Ref. 5b.) Solid square represents the vertex function which for zero external 
momenta gives fmi?- The dotted line in the series for Vm is a dummy line signifying loop closing. 
This line indicates that two different indices are coupled by the dummy momentum p. There are 
two factors in (b). (i) A combinatorial factor 8 from the insertions of A vertices and the subsequent 
noise contraction, and (ii) (2*) for choosing the wave vectors. 
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